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Abstract . This paper summarizes the research on optimal 
trajectories for the National Aerospace Plane performed by the 
Aero-Astronautics Group of Rice University during the period from 
June 22, 1989 to December 31, 1990. It is assumed that the 
aerospace plane is controlled via the angle of attack and the 
power setting. The time history of the controls is optimized 
simultaneously with the switch times from one powerplant to 
another and the final time. The intent is to arrive at NASP 
guidance trajectories exhibiting many of the desirable 
characteristics of NASP optimal trajectories. 

Key Words . Flight mechanics, hypervelocity flight, 
atmospheric flight, optimal trajectories, aerospace plane, 
sequential gradient-restoration algorithm. 
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Notations 

2 

a = acceleration, ft/sec ; 

g e = sea-level acceleration of gravity, ft/sec 2 ; 
h = altitude, ft; 

I S p = specific impulse, sec; 

m = mass,lbf sec 2 /ft; 

M = Mach number; 

q = dynamic pressure, lbf/ft ,* 

Q = heating rate, BTU/ft 2 sec; 

5 = reference surface area, ft 2 ; 

S e = combustor cross-sectional area, ft ^ ; 
t = dimensionless time; 

T = thrust, lbf; 

V = velocity, ft/sec; 

W = mg e = sea-level weight, lbf; 

x = distance along the Earth surface, ft; 
a = angle of attack, rad; 

3 = power setting; 

y = path inclination, rad; 

6 = running time, sec; 

t = final time, sec. 

Subscripts (EMI + EM2) 

0 = beginning of ramjet phase/initial point; 

1 = beginning of scramjet phase; 

2 = end of scramjet phase/final point; 

= final point. 


f 
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Subscripts (EM3) 

0 = beginning of ramjet phase/ initial point; 

1 = beginning of scramjet phase; 

2 = beginning of rocket phase; 

3 = end of rocket phase/final point; 

f = final point. 

Acronyms 

GHAME = general hypersonic aerodynamics model example; 
NASP = national aerospace plane; 

SGRA = sequential gradient-restoration algorithm; 

SSTO = single-stage-to-orbit ; 

TSTO = two-stage-to-orbit . 
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1 . Introduction 

The National Aerospace Plane (HASP) is a hypervelocity 
vehicle which must take-off horizontally, achieve orbital speed, 
and then land horizontally. At this time, its configuration 
not precisely known, but it can be assumed that the powerplant 
includes the combination of four types of engines, turbo j 
turbofan engines for flight at subsonic speeds and low supersonic 
speeds ; ramjet engines for flight at high supersonic speeds; 
scramjet engines for flight at hypersonic speeds; and rocket 

engines for flight at near-orbital speeds. 

in this research, optimal trajectories are studied for a 
given HASP configuration, the so-called general hypersonic 
aerodynamics model example (SHAME) . The optimization study is 
done employing the sequential gradient-restoration algorithm for 


optimal control problems. This algorithm, developed and perfected 
by the Aero-Astronautics Group of Rice University over the years 
1970-85, has proved to be a powerful and reliable tool for 
solving highly constrained/highly nonlinear problems of optimal 
control, such as flight in a windshear and aeroassisted orbital 
transfer. Here, it is applied to the HASP problem. 

in a subsequent research, the current optimization study is 


to be followed by a guidance study. The intent is to develop 
guidance trajectories, capable of approximating the key 
properties of the optimal trajectories. This is to be achieved 
via a feedback control scheme characterized by strong resistance 
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to external disturbances. See, for instance, the guidance schemes 
already developed by the Aero-Astronautics Group of Rice 
University for flight under windshear conditions and aeroassisted 
orbital transfer. 
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2 . Research Results 

The research involves 48 optimization problems obtained by 
combining several performance indexes [PI, P2, P3, P4], 
constraint combinations [A, B, C, D] , and engine models 


[EMI, EM2 , EM3] . 

The performance indexes being minimized are four: 

(PI) weight of fuel consumed; 

(P2) peak dynamic pressure; 

(P3) peak heating rate; 

(P4) peak tangential acceleration. 

The constraint combinations are four: 

(A) Yq = free, q = free, 

(B) Yq — free, q <_ 1500 lbf/ft 2 , 

(C) Yo = 0 *° de< ?' <3 * free, 

(D) Yo = 0,0 de< ?' <3 £ 1500 lbf/ft 2 . 

The engine models are three: 

(EMI) this is a ramjet/ scramjet combination in which the 
scramjet specific impulse tends to a nearly-constant value at 
large Mach numbers; 

(EM2) this is a ramjet/scramjet combination in which the 
scramjet specific impulse decreases monotonically at large Mach 
numbers ; 


a^i = free; 

a T — 2 ^e” 
a>p = free; 

a T < 3g e • 


(EM3) this is a ram jet/ scram jet/rocket combination in 
which, owing to stagnation temperature limitations, the scramjet 
operates only at M £ 15; at higher Mach numbers, the scramjet is 
shut off and the aerospace plane is driven only by the rocket 


engines . 
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Note that a peak heating rate bound, Q 150 BTU/ft sec, is not 
imposed because it can be satisfied or nearly satisfied 
indirectly if the dynamic pressure bound is satisfied. 

The minimization of the above performance indexes is carried 


out under the assumption that the 
completed. The initial conditions 
Xq = 0 ft, 

h 0 = 42000 ft = 12.8 km, 

Vq = 1936 ft/sec, 

Yq = free or Yq = 0*0 deg, 

W Q = 290000 lbf , 

and correspond to Mq = 2, cjq = 10' 
are 

Xf = free, 

h f = 262467 ft = 80.0 km, 

Vf = 25792 ft/sec, 

Yf = 0.0 deg, 

W f = free, 

and correspond to Mf = 27.8, qf = 


turbojet phase has been 
are 

(la) 

(lb) 

(lc) 

(l d) 

(le) 

0 lbf/ft 2 . The final conditions 

(2a) 

(2b) 

(2c) 

(2d) 

(2e) 

11.9 lbf/ft 2 , and orbital 


speed . 

It is assumed that the NASP reference surface area is S = 
6000 ft 2 and that the combustor cross-sectional area of both the 
ramjet and the scramjet is S e = 400 ft 2 . It is further assumed 
that the NASP is controlled via the angle of attack a(t) and the 
power setting B(t), which are subject to the inequalities 
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- 2.0 <_ a £ 12.0 deg, (3a) 

0.0 < 3 < 1 . 0 . ( 3b > 

Finally, it is assumed that the ramjet specific impulse has the 

maximum value I c -r, = 6000 sec for all engine models; that the 

scramjet specific impulse has the maximum value I sp = 2500 sec 

for engine model EMI and I sp = 3100 sec for engine models EM2 and 

EM3; and that the rocket specific impulse has the constant value 

I = 444 sec for engine model EM3 . The rocket maximum thrust is 

sp 

T = 189200 lbf for engine model EM3 . 

For detailed data and results, see Refs. 1-2. A cross 
section of the results obtained is shown in Tables 1-3, each 
containing a different group of problems [Gl, G2 , G3]. In 
analyzing the results of Tables 1-3, the criteria for judging the 
engineering usefulness of the solutions are as follows: 

(Cl) the initial path inclination y 0 should be small to 
avoid overburdening the turbojet engines during the low-altitude 
portion of the flight; 

(C2) the dynamic pressure q should be kept below 1500 
lbf/ft 2 ; 

(C3 ) the heating rate Q should be maintained below 150 
BTU/ f t 2 sec ; 

( C 4 ) the tangential acceleration a^i should not exceed 3g e . 

Group Gl . This group of problems is concerned with the 
effect of the performance index and includes Problems P1A, P2A , 
P3A, P4A. Each of the performance indexes P1-P4 is minimized for 
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constraints of type A and engine model EMI. Clearly, 
unconstrained solutions are obtained, since the initial path 
inclination, the dynamic pressure, and the tangential 
acceleration are free. See Table 1. 

Among the solutions of group G1 , solutions P1A, P2A, P3A 
lead to excessive violation of criteria Cl and C4 , while solution 
P4A is unacceptable in the light of criteria C1-C3. A common 
characteristic of the solutions of group G1 is the steepness of 
the trajectory at the initial point: values of Yg ranging from 
38.3 deg to 50.0 deg are obtained. 

Group G2 . This group of problems is concerned with the 
effect of the constraint combination and includes Problems PlA, 
PlB, PIC, P1D . The performance index PI is minimized for 
constraints of type A, B, C, D and engine model EMI. See 
Table 2. 

Solutions PlA and PlB lead to excessive violation of 
criterion Cl. Unacceptable values of the initial path inclination 
are obtained, specifically, yo = 42.0 deg and yo = 39.4 deg; this 
is not surprising, since Yo I s left free. 

Solutions PIC and P1D are computed for yo = deg; hence, 

they meet automatically criterion Cl. However, solution PIC leads 
to excessive violation of criteria C2-C4. The only solution 
consistent with criteria C1-C4 is solution PlD. 

If one compares the time histories of the dynamic pressure 
and the heating rate for solutions PlA and PlD, one sees that 
they have something in common: a fast initial climb to quickly 
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decrease the air density, so as to contain both the dynamic 
pressure and the heating rate. In solution P1A, this is obtained 
with Yq= 42.0 deg; in solution P1D, this is obtained with a quick 
increase of the path inclination from Yg = 0.0 deg to 
Yg = 25.0 deg. While the unconstrained solution PlA exceeds the 
tangential acceleration bound, the constrained solution P1D 
satisfies the tangential acceleration bound by reducing the power 
setting when the 3g e limit is met. The fuel penalty paid for 
imposing the additional constraints concerning Yg» <3* on ^y 

2% as can be seen by comparing solutions PlA and P1D. 

Group G3 . This group of problems is concerned with the 
effect of the engine model and includes three solutions of 
Problem P1D. The performance index PI is minimized for 
constraints of type D and engine models EMI, EM2 , EM3 . See 
Table 3. 

Clearly, all the solutions of Group G3 satisfy criteria 
C1-C4. In percentage of the aerospace plane weight at the end of 
the turbojet phase, the minimum fuel weight is 34.3% for engine 
model EMI, 44.3% for engine model EM2 , and 60.7% for engine model 
EM3 . Let us assume that the turbojet portion of the flight 
consumes a fuel weight equal to 5.0% of the take-off weight. 

Then, in percentage of the aerospace plane weight at take-off, 
the minimum fuel weight is 37.6% for engine model EMI, 47.1% for 
engine model EM2 , and 62.7% for engine model EM3 . 

The fact that engine model EM 3 carries a severe fuel weight 
penalty is clear. Since the scramjet operation is discontinued at 
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M = 15, the flight portion from M = 15 to M = 27.8 is spent under 

rocket power. Because the largest amount of energy increase takes 

place during the rocket phase and because the rocket engine has 
lower specific impulse , engine model EM3 uses 33% of the fuel 
weight in the ramjet-scramjet phase and 67% in the rocket phase. 
Indeed, the consequences of carrying the oxidizer onboard are 
severe for this type of SSTO vehicle. 

These results indicate that the required fuel weight and 
connectedly the useful payload are dependent heavily on the 
performance of the scramjet powerplant at hypersonic Mach 
numbers. If engine model EM2 is closer to reality, then the SSTO 

mission appears to be feasible. On the other hand, if engine 

model EM 3 is closer to reality, then the SSTO mission appears to 
be marginal at best. 

Of course, improvements are possible in the areas of 
aerodynamic properties and specific impulse properties via 
highly-integrated airframe/engine combinations. Under this 
scenario, the SSTO mission might become feasible. But prudence 
seems to dictate that a TSTO mission deserves concurrent 


consideration . 
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Table 1. Group G1 solutions, effect of the performance index, 
constraints of Type A, engine model EMI. 


Quantity 


Problem 



Units 


P1A 

P2A 

P3A 

P4 A 


(W 0 -Wf)/W 0 

0.337 

0.347 

0.357 

0.550 

- 

max (q) 

1540 

999 

1157 

3751 

lbf/ft 2 

max (Q) 

165 

161 

98 

495 

BTU/f t 2 sec 

max (a T )/g e 

9.1 

5.2 

4.0 

1.1 

- 

^0 

42.0 

50.0 

40.4 

38.3 

deg 

T 1 

34 

54 

48 

144 

sec 

t 2 

409 

475 

731 

704 

sec 

6 f 

443 

529 

779 

848 

sec 

= W 2 and 

e f = 0 2 

for engine model 

EMI . 
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Table 2. Group G2 solutions, effect of the constraint combination, 



minimum fuel 

weight, engine 

model EMl • 



— 

Quantity 


Problem 


Units 



P1A 

PlB 

PIC 

P1D 



(w cr w £ ) /w o 

0.337 

0.340 

0.339 

0.343 

- 


max (q) 

1540 

1112 

1765 

1500 

lbf/ft 2 

— 

max (Q) 

165 

148 

200 

153 

BTU/f t 2 sec 

— 

max(a T )/g e 

9.1 

3.0 

13.7 

3.0 

- 



Y 0 

42.0 

39.4 

0.0 

0.0 

deg 

— 

T 1 

34 

55 

34 

55 

sec 


T 2 

409 

498 

335 

487 

sec 


6 f 

443 

553 

369 

542 

sec 

— 

= V ?2 and 

Cl 

04 

CD 

II 

M-t 

CD 

Dr engine mode! 

L EMI. 




©2 for engine model EMI. 
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Table 3. Group G3 solutions, effect of the engine model, 
minimum fuel weight, constraints of Type D. 


Quantity 


Engine model 


Units 


EMI 

EM2 

EM 3 


<V w f )/w o 

0.343 

0.443 

0.607 

- 

max (q) 

1500 

1425 

1500 

lbf/ft 2 

max (Q) 

153 

157 

110 

2 

BTU/ft sec 

max(a T )/g e 

3.0 

3.0 

3.0 

- 

Y 0 

0.0 

O 

» 

O 

O 

• 

O 

deg 

T 1 

55 

44 

57 

sec 

T 2 

487 

472 

97 

sec 

T 3 

- 

- 

277 

sec 

9 f 

542 

517 

431 

sec 

W f = W 2 and 

0 f = 0 2 for 

engine models EMI, 

EM2. 



Wf = and = 


0^ for engine model EM3 . 
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3 . Abstracts of Publications 

3.1. MIELE, A. LEE, W. Y., and WU, G. D. , Optimal Trajectories 
for an Aerospace Plane, Part 1: Formulation, Results, and 
Analysis , Rice University, Aero-Astronautics Report 
No. 247, 1990. 

Abstract . This report is concerned with the optimization 
of the trajectories of an aerospace plane. This is a 
hypervelocity vehicle capable of achieving orbital speed, while 
taking off horizontally. The vehicle is propelled by four types 
of engines: turbojet engines for flight at subsonic speeds/ low 
supersonic speeds; ramjet engines for flight at moderate 
supersonic speeds/ low hypersonic speeds; scramjet engines for 
flight at hypersonic speeds; and rocket engines for flight at 
near-orbital speeds. 

A single-stage-to-orbit (SSTO) configuration is 
considered, and the transition from low supersonic speeds to 
orbital speeds is studied under the following assumptions: the 
turbojet portion of the trajectory has been completed; the 
aerospace plane is controlled via the angle of attack a(t) and 
the power setting £(t) ; the aerodynamic model is the generic 
hypersonic aerodynamics model example (GHAME) . Concerning the 
engine model, three options are considered: (EMI) this is a 

ramjet/ scramjet combination in which the scramjet specific 
impulse tends to a nearly-constant value at large Mach numbers; 
(EM2) this is a ramjet/ scramjet combination in which the scramjet 
specific impulse decreases monotonically at large Mach 
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numbers; (EM3 ) this is a ramjet/scram jet/rocket combination in 
which, owing to stagnation temperature limitations, the scramjet 
operates only at M < 15; at higher Mach numbers, the scramjet is 
shut off and the aerospace plane is driven only by the rocket 
engines . 

Under the above assumptions, four optimization problems 
are solved using the sequential gradient-restoration algorithm 
for optimal control problems: (Pi) minimization of the weight of 

fuel consumed; (P2) minimization of the peak dynamic pressure; 

(P3) minimization of the peak heating rate; and (P4) minimization 
of the peak tangential acceleration. The above optimization 
studies are carried out for different combinations of 
constraints, specifically: initial path inclination either free 
or given (Yq = 0); dynamic pressure either free or bounded 
(q <_ 1500 lbf/ft^); tangential acceleration either free or bounded 
(a T £ 3g e ) . 

The main conclusions are as follows: 

(a) For an aerospace plane governed by GHAME + EMI, the 
SSTO mission requires a weight of fuel consumed equal to 34.3% of 
the initial weight. 

(b) For an aerospace plane governed by GHAME + EM2 , the 
SSTO mission requires a weight of fuel consumed equal to 44.3% of 
the initial weight. 

(c) For an aerospace plane governed by GHAME + EM3 , the 
SSTO mission requires a weight of fuel consumed equal to 60.7% of 
the initial weight. 
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(d) If one assumes that engine model EM2 is the one 
closer to reality, then the SSTO mission appears to be feasible. 
Obviously, its ability to deliver payloads can be improved via 
progress in the areas of aerodynamic properties and specific 
impulse properties. 

(e) If one assumes that engine model EM3 is the one 
closer to reality, then the SSTO mission appears to be marginal, 
unless substantial progress is achieved in the areas of 
aerodynamic properties and specific impulse properties. Under 
this scenario, alternative consideration should be given to 
studying the feasibility of a two-stage-to-orbit (TSTO) mission. 
3.2. MIELE, A., LEE, W. Y., and WU, G. D. , Optimal Trajectories 

for an Aerospace Plane, Part 2: Data, Tables, and Graphs , 

Rice University, Aero-Astronautics Report No. 248, 1990. 

Abstract . This report is a follow-up to Ref. 1 and 
presents data, tables, and graphs relative to the optimal 
trajectories for an aerospace plane. A single-stage-to-orbit 
(SSTO) configuration is considered, and the transition from low 
supersonic speeds to orbital speeds is studied for a single 
aerodynamic model (GHAME) and three engine models. 

Four optimization problems are solved using the sequential 
gradient-restoration algorithm for optimal control problems: (Pi) 

minimization of the weight of fuel consumed; (P2) minimization of 
the peak dynamic pressure; (P3) minimization of the peak heating 
rate; and (P4) minimization of the peak tangential acceleration. 
The above optimization studies are carried out for different 
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combinations of constraints, specifically: initial path 
inclination either free or given; dynamic pressure either free or 
bounded; tangential acceleration either free or bounded. 
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